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Constraints on possible age spreads within young massive 
clusters in the Large Magellanic Cloud 
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ABSTRACT 

Recent studies have shown that the observed main-sequence turn-off (MSTO) in 
colour- magnitude diagrams of intermediate age (1-2 Gyr) clusters in the LMC are 
broader than would be nominally expected for a simple stellar population. This has 
led to the suggestion that such clusters may host multiple stellar populations, with age 
spreads of 100-500 Myr. However, at intermediate ages, spreads of this magnitude are 
difficult to discern and alternative explanations have been put forward (e.g., stellar 
rotation, interacting binaries). A prediction of the age-spread scenario is that younger 
clusters in the LMC, with similar masses and radii, should also show significant age 
spreads. In younger clusters (i.e., 40-300 Myr) such large age spreads should be readily 
apparent. We present an analysis of the colour-magnitude diagrams of two massive 
young clusters in the LMC (NGC 1856 and NGC 1866) and show that neither have 
such large age spreads, in fact, both are consistent with a single burst of star- formation 
(cr(age) < 35 Myr). This leads us to conclude that either the intermediate age clusters 
in the LMC are somehow special or that the broadened MSTOs are not due to an age 
spread within the clusters. 
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1 INTRODUCTION 

The traditional view of stellar clusters as simple stellar pop- 
ulations (SSP - i.e. all stars have the same age and metallic- 
ity within some small tolerance) has recently been called into 
question due to the discovery of multiple main sequences, 
horizontal branches, turn-offs, and sub-giant branches in the 
old globular (> 10 Gyr) clusters in the Galaxy (e.g., Bedin 
et al. 2004). However, it is currently unclear if this is a gen- 
eral property of massive stellar clusters, or if this feature is 
unique to the globular clusters. High precision photometry 
of intermediate age (1 — 2 Gyr) clusters in the LMC has 
shown that while the main-sequence and giant branches are 
well described by an SSP model, the main sequence turn- 
off is broader than would be expected, even when photo- 
metric errors are taken into account (e.g., Mackey & Broby 
Nielsen 2007, Milone et al. 2009, Goudfrooij et al. 2011a,b). 
This broadening has been taken as evidence for a signif- 
icant age spread within the clusters (> 100 — 500 Myr), 
which if true, could provide vital clues into how the multi- 
ple populations in globular clusters formed (e.g. Conroy & 
Spergel 2011, Keller, Mackey, & Da Costa 2011). 
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Alternatives explanations to age spreads have been put 
forward to explain the extended main sequence turnoff (eM- 
STO) seen in LMC clusters. Yang et al. (2011) have sug- 
gested that interacting binaries can create an eMSTO as 
well as a dual red clump, observed in some intermediate 
age clusters (Girardi et al. 2009). Additionally, Bastian & 
de Mink (2009; see also Yang, Bi, & Meng 2012) suggested 
that stellar rotation may be able to produce an eMSTO due 
to the change in luminosity and colour in rapidly rotating 
stars although this scenario has been called into question 
(Girardi et al. 2011). 

If massive intermediate age clusters do host signifi- 
cant age spreads or multiple populations, this could pro- 
vide a promising avenue to understand the multiple pop- 
ulations observed in globular clusters. Indeed, there have 
been many recent works that have attempted to link the 
suggested age spread in the intermediate age LMC clusters 
with globular clusters (e.g., Conroy & Spergel 2011, Keller 
et al. 2011). These works have suggested that massive clus- 
ters (> lO^M©) clusters have deep enough potentials to re- 
tain material lost due to stellar evolution of the first gen- 
eration of stars and potentially accrete new material from 
their surroundings. If such material can be retained, a sec- 
ond (or multiple) populations of stars can be formed within 
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the existing cluster. This is similar to models put forward 
to explain the multiple populations within globular clusters 
(e.g., de Mink et al 2009, and references therein), although 
such models may require a larger population of first gener- 
ation stars than can be accommodated by observations of 
the number of field stars within the host galaxy (Larsen, 
Strader, k Brodie 2012). 

Other studies have suggested that age spreads may be 
present in young (< 100 Myr), massive (> lO^M©) extra- 
galactic clusters (e.g., Larsen et al. 2011). These potential 
age spreads are significantly less in duration, being tens of 
Myr, instead of hundreds of Myr in the intermediate age 
LMC clusters. However, the authors conclude that interact- 
ing binaries or other effects may explain the results without 
resorting to age spreads, especially as the estimated spread 
in ages appears to be correlated to the age of the cluster 
(Larsen et al. 2011). 

In the present work we look for significant age spreads in 
two massive clusters in the LMC, NGC 1856 and NGC 1866. 
Both clusters are massive (- lO^Mo, see Tablefflfor the ba- 
sic cluster properties) and have properties similar to clusters 
showing the eMSTO feature (see Fig.[T]). If the age spreads in 
the intermediate age LMC clusters are real, we would expect 
to find similar age spreads in these two clusters. Both clus- 
ters have been studied previously, in particular NGC 1866 
has an extensive history of being used for calibration of stel- 
lar evolutionary models (see Rosenfield et al. in prep, for 
an extensive review). Brocato et al. (2003) derived an age 
of - 160 - 250 Myr for NGC 1866 based on an analysis of 
the colour-magnitude diagram (CMD). NGC 1856 has been 
most recently studied by Kerber et al. (2007) who derived 
an age of ^ 280 Myr. We will use the photometry presented 
in each of the above works to search for possible age spreads 
within these clusters. 

Throughout this work we adopt the Padova isochrones 
(Bressan et al. 2012) at a metallicity of z=0.008, a distance 
modulus to the LMC of 18.50 (e.g., Marconi & Clemen- 
tini 2005), and a Salpeter (1955) stellar initial mass func- 
tion. 



2 NGC 1866: HESS DIAGRAMS AND THE 
STAR FORMATION HISTORY 

NGC 1866 is a relatively young (100-200 Myr), massive clus- 
ter (-- lO^Mo) in the LMC (e.g. McLaughlin & van der 
Marel 2005). In order to constrain a possible age spread 
within this cluster we use Hubble Space Telescope - Wide 
Field Planetary Camera 2 (HST-WFPC2) images taken in 
the and F8 l^W h^nds. The data were presented in 

detail in Brocato et al. (2003), who analysed the images to 
study the stellar luminosity function within the cluster in or- 
der to constrain the stellar mass function and the efficiency 
of core overshooting in stellar models. In the present work 
we use their stellar photometric catalogues that have been 
transformed to the standard VI bands. 

Here, we focus on the inner 4.8 pc of the cluster, in or- 
der to avoid background contamination as much as possible. 
No background subtraction was performed, as background 
stars will only increase the appearance of any potential age 
spreads that are present. 

In Figure [2] we show the CMD of the observed 
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Figure 1. The mass vs. effective radius for a sample of rela- 
tively young (< 500 Myr) massive (> lO^M©) clusters in the 
LMC (black dots), taken from the catalogue of McLaughlin & 
van der Marel (2005). The two clusters discussed in the present 
paper are marked with large open circles. Intermediate age clus- 
ters that show the eMSTO feature, from the sample of Goudfrooij 
et al. (2011) are shown as filled (red) triangles. Note that the clus- 
ters presented in the current work are more massive than many of 
the clusters showing the eMSTO feature and have similar effective 
radii. 

stars (back points) along with contours (solid dark lines) 
representing the number density of stars at a given 
colour /magnitude. The main sequence is clearly evident and 
the helium burning stars have two over-densities. Addi- 
tionally, in the six panels we show Hess diagrams adopt- 
ing the Padova isochrones and our photometric errors for 
six different ages, ranging from 100 to 250 Myr. We find 
A\ — 0.15, consistent with previous studies (e.g., Brocato 
et al. 2003). Focussing on the position and morphology of the 
helium burning stars, we see that the best fit is obtained at 
180 Myr. Significant numbers of stars from younger ages 
can be excluded as they would produce helium burning star 
features that do not correspond to the observations. Older 
populations can also be excluded due the lack of stars with 
similar colours, but fainter than, the observed helium burn- 
ing stars. 

One feature that should be noted is that the observed 
main sequence extends to brighter magnitudes than the 
predictions of the Hess diagram for the best fit age, by 
^ 0.5 mag. These stars could be due to a younger gen- 
eration, although evidence against this is seen in the lack 
of corresponding post-main sequence (i.e., helium burning) 
stars. Alternatively, these stars could be due to binaries in 
the cluster and/or crowding in the images since the observa- 
tions are focussed on the relatively dense core of the cluster. 

In order to provide a more quantitative measure of pos- 
sible age spreads within NGC 1866, we also fit the star- 
formation history (SFH) of the cluster using the observed 
CMD, Padova isochrones and the FIT SFH code presented 
in Silva- Villa & Larsen (2010) and Larsen et al. (2011). We 
restricted the fits to brighter than My = 1 in order to avoid 
incompleteness issues. We refer to the reader to the above 
works for details on the fitting procedure. Photometric errors 
were included in the fits, and we chose boxes in the CMDs 
(regions within the CMDs to compare the observations with 



© 2013 RAS, MNRAS 000, [T]-?? 



A lack of age spreads in young LMC clusters 3 

ID log (Age/yr) log (Mass/M©) Reff (pc) Rcore (pc) Vesc (km/s) 

NGC 1856 8.12 (8.45) 4.88 7.8 1.7 10.5 

NGC 1866 8.12 (8.25) 4.91 10.4 2.8 9.2 



Table 1. Properties of the two LMC YMCs discussed in the present work. All values are taken from McLaughlin & van der Marel (2005), 
adopting King (1966) profile fits. The number in parenthesis is the best fit age from the star formation history analysis presented in this 
work. Note that Vesc would be a factor ~ 1.2 — 1.5 higher if we applied the corrections of Goudfrooij et al. (2011b) to estimate Vesc at 
an age of 10 Myr. 




Figure 2. Hess diagrams for NGC 1866. The dots represent individual stars, the black contours show the binned data (i.e. the number 
of stars in a given colour/magnitude bin) and the filled (colour) contours represent the models at the age given in each panel. Note the 
strong constraints given by the morphology /posit ion of the helium burning stars. The best fitting age is ~ 180 Myr, with an acceptable 
range of ~ 30 Myr on either side, based on the helium burning stars. Note the handful of stars above the nominal main sequence for the 
best fit age. 



the isochrones) with the following limits: for NGC 1866 
-0.3 ^ {V - I) ^ 0.3 and -4.0 ^ My ^ 1.5 for the main 
sequence and 0.3 ^ {V - I) ^ 1.5 and -4.0 ^ My ^ 0.05 
for the post-main sequence. We experimented with differ- 
ent limits and found that the results were not sensitive to 
the exact choice of colours/magnitudes over which to fit. 
The resulting SFH is shown in Fig. O under the assump- 
tion that no binary stars are present, as solid circles. The 
solid (blue) lines show two Gaussians fit to the data. The 
main peak has a dispersion of 18.4 Myr. The secondary peak 
is due to the handful of stars above the nominal main se- 
quence turnoff (for the best fitting age). This peak can be 
reduced or removed by including the effects of binaries in 
the fitting procedure, although crowding effects are likely to 
be contributing as well. However, the inclusion of binaries 
increases the dispersion of the main peak from 18 Myr to 
30 Myr. 



From this we can see that any possible age spread is 
limited to 30 Myr of the peak age. This is significantly less 
than the 100-500 Myr expected from observations of the 
intermediate age clusters in the LMC. 



3 NGC 1856: HESS DIAGRAMS AND THE 
STAR FORMATION HISTORY 

Another of the most massive young (< 1 Gyr) clusters in the 
LMC is NGC 1856. Using the mass and size estimated in the 
survey of Macky & Gilmore (2003), Goudfrooij et al. (2011b) 
suggest that this is the only cluster with an escape velocity, 
Vesc, larger than the empirically derived "critical" escape 
velocity of 10 km/s, hence that it should host multiple gen- 
erations of stars within it. 

We can test this suggestion by studying the CMD of 
NGC 1856 presented in Brocato et al. (2001). These au- 
thors present HST-WFPC2 F450W and F555W observa- 
tions of the cluster, with the PC chip positioned near the 
cluster core. We refer the interested reader to that paper 
for the observational details. The CMD, uncleaned for field 
and foreground star contamination, is shown in Fig. [4l with 
points representing individual stars and the solid (black) 
line contours show the number density of points at that 
colour /magnitude. As in §[2] we restrict our analysis to the 
central 4.8 pc to minimise contamination of the sample. We 
also show theoretical Hess diagrams for three different ages 
in the panels. 

Using the position and morphology of the helium burn- 
ing stars we can constrain the age and any age spread present 
within the cluster. From Fig. Owe see that the helium burn- 
ing stars are best reproduced by the 280 Myr isochrone. 
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Figure 3. The derived SFH of NGC 1866 (solid cirlces) and 
NGC 1856 (open circles), assuming no binary stars. The points 
represent the estimations, while the solid (blue) and dashed 
(red) lines show Gaussian fits to the points for NGC 1866 and 
NGC 1856, respectively. For NGC 1866 the main peak has a 
dispersion of cr = 18.4 Myr and a peak at 180 Myr, while for 
NGC 1856 we find a dispersion of 33.6 Myr and a peak at 280 Myr. 
The dispersions are upper limits to the true age distribution. The 
secondary peaks are driven by the handful of stars above the main 
sequence turn-off (see Figs. [2]&;[4|), which are likely due to bina- 
ries and crowding effects in the central regions of the clusters, 
since no corresponding post-main sequence features are seen. 

However, due to the lower signal-to-noise of the data, the 
constraint on the possible age spread within the cluster 
(based on helium burning star morphology) is weaker than 
in the case of NGC 1866. We can rule out significant pop- 
ulations of stars with ages less/more than 50 Myr than the 
best fit isochrone. 

As in NGC 1866, there are a handful of stars above the 
nominal main sequence turnoff. These are presumably bi- 
nary stars and/or stars affected by crowding due to the high 
surface density of stars in the central parts of the cluster. 
The possibility of these stars representing a younger pop- 
ulation is low, due to the lack of corresponding post-main 
sequence stars. 

We have also estimated the SFH of NGC 1856 using the 
same method as in § [2] For this cluster we adopted fitting 
boxes with limits of -0.3 ^ (B - V) ^ 0.3 and -4.0 ^ 
Mv ^ 0.0 for the main sequence and 0.3 ^ {V — I) ^ 1.5 
and —4.0 ^ My ^ 0.0 for the post-main sequence. We find 
Ay — 0.8, consistent with previous studies (e.g., Kerber 
et al. 2007), and the estimated SFH is shown in Fig. [3] as 
open circles. The dashed lines show two Gaussian fits to the 
estimated distribution. We find a best fit age of 280 Myr 
and a dispersion of ^ 35 Myr. As in the case of NGC 1866, 
the secondary peak at younger ages is due to the handful of 
stars above the nominal main sequence of the best fit age, 
presumably due to binaries and crowding effects. 



4 DISCUSSION AND CONCLUSIONS 

We have tested the assertion that large age spreads (> 
100 Myr) are a common feature of massive clusters. These 
suggestions were based on previous studies of massive (> 



lO^M©), intermediate age (1-3 Gyr) clusters that showed 
that the observed main sequence turnoff (MSTO) was 
broader than predicted by simple stellar population, which 
could be explained by an extended star-formation history. 
Alternatively, it has been suggested that extended MSTOs 
(eMSTO) could be due to stellar rotation or interacting 
binaries. If the eMSTOs are due to an age spread within 
the cluster, then it would be expected that younger clusters 
(ages < 500 Myr) with similar properties (mass and radius) 
should show a clear age spread (100 — 500 Myr) within them. 
Such an age spread would be easily detectable through the 
analysis of their CMDs. We note that the eMSTO clusters 
are generally better described by a continuous distribution 
of ages, rather than a bimodal distribution (e.g., Goudfrooij 
et al. 2011b). 

Using CMDs based on HST imaging of two massive 
young clusters in the LMC (NGC 1856 and NGC 1866; see 
Table [T] for their properties) we find that both are consis- 
tent with a single burst of star formation, and place upper 
limits to possible age spreads within them of cr = 18 Myr 
and cr = 35 Myr for NGC 1866 and NGC 1856, respectively 
(fitting the derived SFH with a Gaussian). These are upper 
limits to the true dispersion, as crowding, contamination, 
and differential extinction has not been taken into account 
(additionally the errors in the observations introduce some 
minimum width to the best fit SFH). The lack of an age 
spread in each of these clusters is consistent with that found 
for massive (10^ — lO^M©) young (< 10 Myr) clusters in the 
Galaxy, such as NGC 3603 and Westerlund 1 (e.g. Kudryavt- 
seva et al. 2012). In Fig. [5] we show Gaussian fits to the SFH 
of both clusters, as well as the expected distribution based 
on the inferred SFH of seven intermediate age clusters pre- 
sented in Goudfrooij et al. (2011b), which have been shifted 
to an age of 200 Myr for comparison. The two clusters pre- 
sented in the current work clearly do not match the inferred 
SFH of the intermediate age clusters. 

The lack of age spreads within the two young clusters in 
the current work can be interpreted in two ways. Either the 
intermediate age clusters with eMSTOs are fundamentally 
different than the clusters presented here in some intrinsic 
property (although their masses and radii are comparable) 
or the eMSTO feature observed in these clusters is not due 
to an age spread. Alternative explanations for the eMSTOs 
have been put forward in the literature, namely interacting 
binary stars (Yang et al. 2011) and stellar rotation (Bas- 
tian & de Mink 2009). Girardi et al. (2011) have called the 
rotation scenario into question, by correctly pointing out 
that the models presented by Bastian & de Mink were not 
isochrones, as the influence of rotation on stellar lifetimes 
was not taken into account. When Girardi et al. did cor- 
rect for the changing lifetimes, the spread in the modelled 
MSTO decreased. However, we note that the rotation rate 
used by Girardi et al. was rather low (using the mean of the 
Royer et al. (2007) rotational distribution and not taking 
into account the significant numbers of more rapidly rotat- 
ing stars) and that their models did not show any difference 
between the rotating and non-rotating stars on the main se- 
quence, although rotation should affect stellar structure and 
be observable on the main sequence (e.g., Brott et al. 2011; 
Ekstrom et al. 2012). 

Both of the clusters discussed here are well above the 
lO^M© limit suggested by Conroy and Spergel (2011) in or- 
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Figure 4. Similar to Fig.[2]but now for NGC 1856. Using the helium burning stars position and morphology, the best fit age is ~ 280 Myr. 
Due to the noisier data, the age spread constraint is less than in NGC 1866, with significant populations of stars with ages ±50 Myr not 
consistent with the observations. As in NGC 1866, there are a handful of stars brighter than the nominal main sequence turnoff, for the 
best fit age, presumably due to crowding and binaries, as no counterpart post-main sequence population is seen. 
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Figure 5. A comparison between the estimated SFH of 
NGC 1856 and NGC 1866 in the current work (solid blue lines) 
and the expectations from seven intermediate age clusters in 
Goudfrooij et al. (2011b) shifted to a reference age of 200 Myr. 
All clusters have been normalised so that the peak contribution 
is unity. The distributions of the two clusters presented in the 
current work are significantly narrower, hence do not conform to 
the expectation if age spreads were a common feature in massive 
clusters. 



der for clusters to host multiple populations. Both clusters 
have an escape velocity of > 10 km/s (after corrections for 
mass loss and change of radius, using the method discussed 
in Goudfrooij et al. 2011b) which is larger than many of 
the clusters that display the eMSTO feature, and the clus- 
ters discussed here have similar properties as the eMSTO 
clusters. 

The observations presented here, of a lack of an age 
spread within two young massive clusters, suggests that the 
eMSTO feature in intermediate age clusters is not due to age 
spreads. This is also consistent with the lack of abundance 
spreads within intermediate age clusters (e.g., Mucciarelli 
et al. 2007), and suggests that the eMSTO feature is not 
related to the multiple stellar populations seen in globu- 
lar clusters. Additionally, the observed colour spread in the 
eMSTO feature is similar in all intermediate age clusters 
observed to date, resulting in an age spread that is propor- 
tional to the age of the cluster, suggesting that other effects 
may be causing the eMSTO feature. 



In a future work, we will present an analysis of the ob- 
served CMDs of another ^ 15 young LMC/SMC clusters. 
Additionally, we will present modelling of each of the clus- 
ters, including the two clusters presented here, that involves 
a more realistic treatment of the errors as well as the effects 
of crowding and binarity. Through this method we will be 
able to place the tightest possible constraints on possible 
age spreads within these clusters. 
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